Current behavior of a quantum Hamiltonian ratchet in resonance 
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We investigate the ratchet current that appears in a kicked Hamiltonian system when the period 
of the kicks corresponds to the regime of quantum resonance. In the classical analogue, a spatial- 
temporal symmetry should be broken to obtain a net directed current. It was recently discovered 
that in quantum resonance the temporal symmetry can be kept, and we prove that breaking the 
spatial symmetry is a necessary condition to find this effect. Moreover, we show numerically and 
analytically how the direction of the motion is dramatically influenced by the strength of the kicking 
potential and the value of the period. By increasing the strength of the interaction this direction 
changes periodically, providing us with a non-expected source of current reversals in this quantum 
model. These reversals depend on the kicking period also, though this behavior is theoretically 
more difficult to analyze. Finally, we generalize the discussion to the case of a non-uniform initial 
condition. 
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I. INTRODUCTION 

The ratchet effect, a directed transport without any ex- 
ternal net force, has attracted an increasing interest since 
early studies by Feynman et al.[l|. This phenomenon 
has a wide range of possible applications in rectifiers, 
pumps, particle separation devices, molecular switches, 
and transistors (see review articles [2|, |3[ and the refer- 
ences therein) . It is also of great interest in biology, since 
the working principles of molecular motors can be con- 
veniently explained in terms of ratchet mechanisms 
Moreover, it is possible to demonstrate quantum ratchet 
effects by using cold atoms. As a result, many dif- 

ferent scenarios have been considered [§. [9l ITol. [TH Ha [l3{ ■ 

At the classical level, directed transport in periodic 
systems can be associated to a broken spatial-temporal 
symmetry [l4]. We can obtain a net current by means 
of a periodically kicked system, for instance. In this sit- 
uation one can break the spatial symmetry by using an 
asymmetric potential and the temporal symmetry by in- 
troducing^ dissipation [IH or an asymmetric kicking se- 
quence [lOf . In these cases the quantum versions present 
the same symmetry features as the classical counterparts, 
showing the corresponding current. 

Current reversals is one of the interesting ratchet fea- 
tures that has attracted considerable interest [n| [T3, [H, 
[H, H3] • Besides the ones due to symmetries in the po- 
tential, there are essentially two types of current rever- 
sals. One corresponds to dissipative chaotic systems and 
it is originated by bifurcations from chaotic to periodic 
regimes, that is, by the transition from strange attractors 
to simple ones 



17] . The other [20( is explained by the 



rected transport associated to spatial asymmetry plus 
quantum resonance effects rather than to an explicit 
spatio-temporal symmetry breaking. This kind of sys- 
tems was introduced in [2l| where the authors have found 
a new mechanism for directed motion in quantum Hamil- 
tonian systems. In this kind of systems momentum grows 
indefinitely, i.e. it does not stabilize around an asymp- 
totic value (we could see this as a "rectification of force" ). 
In their work they have shown that even if the system is 
time symmetric there can be transport in quantum res- 
onance. Quantum resonance (QR) is a pure quantum 
phenomenon without a classical counterpart. In the well- 
known kicked rotor (KR) system, 
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fact that below certain temperatures, quantum tunneling 
can cause a change in the direction of transport. 

The system that we consider in this paper shows di- 



where St = ^2 n S{t — nT), T is the period of the kick, 
k the strength of the kick, and H, the Planck constant, 
and the moment of inertia of the rotator has been taken 
equal to 1). Given a value of the kick strength k, spe- 
cial resonant regimes of motion appear for periods with 
values T = 4-7r^, where the integers r and q are mutu- 
ally prime. Under these conditions the system regularly 
accumulates energy which grows quadratically with the 
time and with k [22|], that is (p 2 ) cx k 2 t 2 , where p stands 
for the momentum operator (for simplicity we refer to 
(p) as the momentum henceforth). 

In the quantum KR at resonance (and for a symmetric 
or antisymmetric initial condition) there is no growth of 
the momentum, i.e., it is always equal to its initial value 
(though it does grow like (p) cx kt if the initial condition 
is generic). In order to find a net current for any initial 
condition, a different potential has to be used. We prove 
that breaking any spatial symmetry is a necessary con- 
dition. The potential that we study corresponds to the 
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double well-kicked rotor (dw-KR) [10| . This potential has 
been experimentally realized in optical lattices (23j . 

In this model it is possible, for high resonances (q > 2), 
to have directed transport even for symmetric initial 
conditions. In addition, this modified KR at resonance 
shows a new kind of current reversals that (although be- 
ing of quantum origin) are not due to tunneling (like in 
[2fJ). In fact, in our model the momentum will evolve 
as (p) oc g(k)t where g(k) is a non monotonic function 
of k. This is quite different from that case in the usual 
KR model. Moreover, we are able to give analytical pre- 
dictions for the current reversals following a perturbative 
approach. 

The paper is organized as follows. The main part of 
our paper, Section II, is devoted to the numerical and 
analytical study of the behavior of the current, including 
current reversals due to the variation of the kick strength 
and period. We study this phenomenon analytically for 
small values of the asymmetry. We also generalize this 
discussion to non-uniform initial conditions. In Section 
III we show that breaking the spatial symmetry is a nec- 
essary condition to find directed current at quantum res- 
onance. In section IV we present our conclusions. 



II. DIRECTED CURRENT BEHAVIOR 



The Hamiltonian of the system is given by: 
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H = — 2W2 +k{cos{6) 



acos(2# + a))8i 
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where St = J2 n " T 1S * ne P er i°d and k is the 
strength of the kick. This Hamiltonian can be used to 
study a gas of cold atoms in an optical lattice. The dif- 
ference between the dw-KR and the KR is due to the 
parameters a, that is the relative strength of the second 
harmonic and a, which is a parameter that breaks the 
symmetry of the kicking potential (V(6) ^ V(2tt — 6) 
where V{9) = fc(cos(0) + acos(26> + a))). The QR con- 
dition is given by T = 4irr/q, as in the KR case. Being 
h = 1, an effective Planck's constant can be defined as 
fi e ff = 8w^T = T. It is important to underline that 
there is no dissipation in our model. We study the prob- 
lem on a torus and the initial condition is 4>o — 1/ \/2tt, 
unless otherwise stated. This can be implemented on 
an optical lattice with a wavefunction almost uniformly 
spread on many lattice sites. The choice of this initial 
condition is due to the fact that it is symmetric. This 
leads to the non typical behavior of this system compared 
with the usual KR where for this scenario the momen- 
tum would stay always at zero. Incidentally, it makes the 
analytical computations much easier. 

Firstly, we would like to make reference to the inset of 
FigJTJ where we can see that for a given fc, the momentum 
(p) grows linearly with time, that is with the number of 
kicks N ((p) oc t oc N). For this numerical evaluation 
we chose T = An 1/3, k — 5, a = 0.01 and a = tt/4. 
Since (p) grows linearly with the number of kicks N we 
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FIG. 1: The average effective force < / > as a function of the 
kick strength k, for r/q — 1/3, a = 2 and a = n/4. In the 
inset we show the momentum (p) versus the number of kicks 
for k — 5, r/q = 1/3, a = 0.01 and a = tt/4. 



focus our attention on a quantity defined as (/) = (p) /N, 
which we call average effective force. In FigQ]we show 
the average effective force versus the kick strength. We 
can see that there is a general growth of (/) with fc, but 
more interestingly, we also see that (/) oscillates radi- 
cally, going from positive to negative values, i.e., there is 
current inversion. To obtain this set of values we have 
used a = 2, representative of the behavior of (/) for large 
values of a, and a = tt/4. In order to see if this behavior 
can be detected in an experiment with cold atoms, we 
have evaluated (p) / -J (p 2 ) — (p) 2 . For this set of param- 
eters, we have found that it saturates to the value 0.18 
at about 15 kicks (N = 15), guaranteeing the feasibility 
of a future realization. 

In the following we will show the origin of the current 
inversions that occur for different values of the kick inten- 
sity or the period. To do this we do a perturbative study 
for small values of the parameter a, which is the ampli- 
tude of the second harmonic. In Ref. 22| it has been 
shown that the one kick evolution of an initial condition 
cj> (6) is given by: 
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<M0) = Z>n<M0 + — ) 

n=0 ^ 
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where So n are given by /?o(#)7n with (5 n = exp(— iV(9 
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with a uniform initial condition we have: 

M0) = Po Ef =o + 2 ~it) = ik 

_ ' (4) 

So, after a generic kick N the wavefunction 4>n+i is given 
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by: 



= PoJ21 u ...l N =0^hP(h+h + ...lN) X Tl2p[l 2 +...l N )— X IInPin 
= A) Ea Ilm,n=l 7i m A„ 

.(5) 

where A stands for the combinations given in the previous 
line of this equation. 

The momentum at time N + 1 is given by (pn+i) = 
Jo^W+iHg^+i)^ (ft = 1) giving: 

(PN+1) = 5* (Pq T,A Ilm,n=l 7* m A*J rfi A 

which can be rewritten as: 

(Pn+i) = a* iT (ft* Ea n^,„=i 7f m A*J 
Ea IlZ,n=r Tin, A„ + A) Ea (i ( 7i m 

(7) 

As we have shown in the inset of FigJT] the momentum 
grows linearly with N. The only part of the previous 
expression that satisfies this condition is: 

(PN+l) « £ iT (&* Ea d 7LC) ^ 
(AEi (l(llm,«=lTiA)))^ 

because we have the derivative of a product of N terms. 
Hence, we can approximate Eq.([S]) by: 



(p N+1 ) « ( ^ ( E m# „ ^PM' n )d9 + 

+ ST {Y,l, m ,n,slll* m lnlsPi +m P* m {l3 n+ sl3 s )')de + 



L m . n = ka Nj* n <y n 
( sin(^n -w m>n )x 

ios + a - 2a; ro , n ) - cos (*E2 + a _ 2u; ro , n ) 

x(fe Ji(w m , n A;) - 2 J 2 K,n*:)/um,»)+ 
-cos (^n _ Wm ,„)x 

3 in (4sa + a _ 2w m , n ) - sin + a - 2w m ,„) 

X2 J2(u m ,nk)/u m , n + 

a - 2uj m . n )j 2 {£l m ,nk) 



-2sin(^- 
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where fL„ „ = 



Here we have used H m ,n = [cos(27rn/q) — cos{2irm/q)] 
and f m ,n = [sin(27rn/g) — sin(27rm/q)]. The momentum 
(p) depends on the period parameter r through the coef- 
ficients 7 m . We now see that it is not surprising to find 
current inversion. In fact, we expect a different sign of 
the average effective force for different values of k, be- 
cause the Bessel functions oscillate. 

From the properties of the Bessel functions we can un- 
derstand the behavior for small and large k. For the 
case of k 3> 1, we must guarantee that ka <C 1. In 
general we have that for x <C 1, J a {x) ~ r(i+a) (f ) 
(r is the function which generalizes the factorial to non- 
integer numbers) and for x S> 1 we have that J a (x) ~ 



—J cos(a;— q/4— 7r/2). For higher values of the period 

parameter q (T = An r/q), the situation is more complex 
because we have many different Q m ,n and so there is a su- 
perposition of many terms like cos(r2 m , n fc — air/4 — tt/2). 
This is much more difficult to study analytically. 

In general, for large values of k, we can approximate 
the behavior by 



(9) 

Here we make a further approximation, we drop all the 
integrals with the exception of the first one. This is jus- 
tified by the fact that their integrands are highly oscil- 
lating functions of 9 compared to the first one, so their 
contribution to the final result is negligible. 

In order to simplify the analytic treatment without los- 
ing the essential features of (p) that we want to describe, 
we focus on the usual KR perturbed by a second har- 
monic, i.e., we take k ~ 0(1) and a <C 1 . In this case, 
as long as ka -C 1 we can approximate (3 m ~ (3 m = 
exp(— ikcos{9 + 2irm/q)){l + ikacos(29 + Anm/q + a)). 
Then, Eq.© becomes: 
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(10) 

These terms are all integrable and after a few com- 
putations we can write the momentum as (pn+i) ~ 
E m ,„ L m.n where L m . n is given by: 
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COS(Q m ,„A; - ) + 



where A and B are two constants that can be obtained 
from Eq. (fTT|) . We can think of this expression as divided 
into two parts. The first part dominates for large val- 
ues of k when it is far enough from its zeros. Even if 
we consider only the first part, the sum of many cosine 
functions with different periods gives a very fluctuating 
behavior. 

To confirm the above analysis, we compute the partic- 
ular cases r/q = 1/3 and r/q = 1/5 and compare them 
with the numerical solution obtained from the evolution 
of the wavefunction. In the case of r/q = 1/3, Q m ,n takes 
only two values, either or y3 and so we only have terms 
like J^s/Hk) where rj = 1, 2. 

After some analytical computations we can see that 
the average effective force is given by: 



(/) = k a sin(a) 
- l) k J l (V3k) + | (1 + V3) J 2 (V3fe) 



(13) 
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FIG. 2: Average effective force {/} versus k for a = 0.01 
and a — n/3. Results from the numerical evolution of the 
wavefunction (circles) and the analytical approximation (solid 
line) are compared. In (a) T = 4ir 1/3 while in (b) T — 
47r 1/5. We can see the oscillations showing current reversals. 



The case for r/q = 1/3 is shown in Fig[^a) where the 
numerical (circles) and analytical (solid line) results are 
compared and show good agreement. A different period 
implies different values of fijj and so different values of 
the period of oscillation for the same value of k. We 
should expect very different behavior for different r/q. 
This is clarified in Fig[^b) where we plot the numerical 
(circles) and analytical (solid line) results showing the 
average effective force versus k for r/q = 1/5. In this 
case, differences are due to the cut-off of oscillating terms. 

It is interesting to notice one important difference in 
this model compared to the usual quantum KR model. 
In our model the average effective force does not grow 
linearly with the intensity of the kick. Instead, it oscil- 
lates periodically. Even the peaks do not show this linear 
behavior with respect to k. There are two components 
contributing to the average effective force. One of them 
causes the peak to grow with a power 1/2 ((/) oc fc 1 / 2 , 
i.e. (p) cx A; 1 / 2 *) and the other makes it to grow with a 
power 3/2 ((/) cx fc 3 / 2 ). 

Reversals due to changes in the value of the kick 
strength k have been found previously [24| , but in a com- 
pletely different context. Here we find reversals due to 
changes in k and also in the value of the period, i.e. in 
the effective Planck constant. Since in our case the na- 
ture of the current is different from the one in [24j , these 
inversions are also of a different and novel character [l3| . 
An insightful model of QRs in the quantum KR can be 
found in (25| . and we think that a theoretical explana- 
tion for this phenomenon could be investigated following 
those lines. 

In the generic case, the initial condition will not be 
symmetric, a situation which is most likely to happen in 
real experimental conditions. As already discussed, for 



an asymmetric initial condition and a symmetric kick, 
quantum resonance induces directed motion and the mo- 
mentum will increase linearly with the time and with the 
strength of the kick. On the other hand, with a non 
symmetric kick in quantum resonance, we will see the 
oscillating behavior that we have found previously, su- 
perimposed in this case to a linear growth due to the 
asymmetry of the initial condition. This is confirmed by 
the numerical results shown in Fig. [3l We have con- 
sidered T = 4ir 1/3, a = 0.01, a = tt/3 and the initial 
condition (j>o(8) — ?7COs(cos(#) + sin(26>)), where r] is the 
normalization constant. 
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FIG. 3: Average effective force (/) versus k for T = 4n 1/3, 
a = 0.01 and a = 7r/3. The initial condition is <j>o{ff) = 
77Cos(cos(#) + sin(20)), where r\ is the normalization constant. 



III. NECESSARY CONDITION FOR DIRECTED 
TRANSPORT IN QUANTUM RESONANCE 

In this section, we shall show that in the absence of 
a net force, either an asymmetry in the initial condition 
or in the kicking potential is necessary to have directed 
transport. It should be stressed that these are neces- 
sary but not sufficient conditions, because it is possible 
to have zero transport also when at least one of these 
two asymmetries is present. Representing a gas of cold 
atoms with the wavefunction tf>, the momentum is given 

by (p) — —ifi ^* iwfyd®- Let's take again a generic 
kicked system with a Hamiltonian given by: 

H = -\w +m&r (14) 

where St = Tl n S{t — nT), T is the period of the kick 
and V is the shape of the kicking potential. We assume 
that the system is at quantum resonance (T — At:- with 
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r and q relative prime numbers) . Given the wavefunction 
after N kicks (</>n), we can obtain the wavefunction after 
N + l kicks 4> N + 1 by Q: 



(15) 



exp(-iV(0 



where /3 n 

Erio ex p(-^ 

7 one can verify that 

In = In+q = lq-n- 



-)) and j v 



-). Using the definition of 



(16) 



Let's study the case with a symmetric initial condition 
(4>n(2it — 8) = 4>n{9)) and a symmetric kicking potential 
(V(2ir — 9) = V(6)). In this case we can see that /3„(27r — 
9) = (3 q - n (9) and then 

<j) N+1 (2Tr -6) = /3 (2^ - 9) J2 n ln cp N (2n - 9 + 2n^ 
n+1 (2tt -6)= po(0) E„ J q -n<pN (0 + 2vr^Ii) . 

This shows that given a symmetric initial condition and 
a symmetric kicking potential, the wavefunction will al- 
ways be symmetric and therefore the momentum will al- 
ways be zero. This result is readily extended to anti- 
symmetric initial conditions. 

We now show that if the initial condition has a shape 
like <t> N {9) = e iLe f N {6) where f N {9) = f N {2n - 9) is 
even, the kick will not change the value of the momentum 
which is (p) = L. For example, an eigenstate of the 
momentum O = ^7= wm n °t change its direction or 



speed. Using Eq. (fT5|) it is easy to show that: 
<Pn + i{9) = 0o(O) Enlne lL{e+2 ^ ) f N (9 + 2n\ 



= M9y Le Y. n ^ 



N 



2tt^ 



. (18) 

where we have factorized out the e lLe which contributes 
to the momentum and now the factors 7„ are multiplied 
by a phase e lLI > 27T ~v\ The symmetry of 4>n+i(6) is given 
by: 



N+1 (2n-9) = 
A)(2tt - 9) £„ Jn e iL ^- e+ ^V f N (2tt 



9 + 2-k'- 



Po(6)e 



iL(2-!r-9) 



iL(27rS) 



f N (9 + 2n^ 



/3 o (0)e^-*)£ n7g -ne iL(2 ^ J /. 
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(19) 



Using the fact that 7„e 



q = e i'j n+ 2L = 

^ 2 (2-f K . „ = ^ L2 ^ lq ^ 2L it is possible to 



e q 'l q -n+2L = e 

show that: 
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iL(2ir^ 
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(20) 



From Eq. (|20|) it is evident that after any kick 4>n(0) can 
be decomposed as 0at(0) = e lLe /jv(0) and this proves 
that the momentum will not change. 

It is well known that with an asymmetric initial con- 
dition we can have transport. In fact, for the simple case 
T = 4tt, we have that <p N = e~ lArv ^Vo- If <j>o is not 
symmetric, even if we have zero net force the momentum 
grows linearly with the number of kicks. In [2lj . it has 
been found that with a non-symmetric kick, transport 
can be induced by high q resonances. This means that 
it is necessary to at least break the spatial symmetry to 
have transport. To do this we can either start with a non- 
symmetric initial condition, use a non-symmetric kick or 
a combination of the above situations. 



IV. CONCLUSIONS 

We have analyzed in detail the phenomenon of directed 
transport for a Hamiltonian system in quantum reso- 
nance. We have found that directed transport can be 
in different directions depending on the intensity of the 
kick or on its period. This phenomenon was unexpected 
and presumably due to how the intensity of the kick and 
the value of the period affect the gradient of the phase 
of the wavefunction. Very interestingly, the direction of 
the motion changes periodically with the intensity of the 
kick. Also remarkable is the fact that the momentum at 
the peaks does not increase linearly with the intensity of 
the kick, but it follows one term that goes as a square 
root and another that goes as a power law with exponent 
3/2. We have also shown that for an asymmetric initial 
condition the periodic effect due to quantum resonance is 
superimposed to the expected drift due to the asymmetry 
of the initial condition. Finally, we have found that even 
though an explicit time symmetry breaking is not needed, 
breaking the spatial symmetry is a necessary condition 
to have directed transport. By numerical simulations we 
have also found out that this effect is very sensitive on 
the initial condition. To see clearly these effects exper- 
imentally a non-interacting BEC could be used because 
the initial wavefunction is less-spread in the momentum 
space than for a gas of cold atoms. 
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